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Abstract 

There is event by event geometric as well as quantum fluctuations in the initial condition of heavy-ion collisions. The standard 
technique of analysing heavy-ion collisions in bins of centrality obtained from final state multiplicity averages out the various 
initial configurations and thus restricts the study to only a limited range of initial conditions. In this paper, we propose an 
additional binning in terms of total spectator neutrons in an event. This offers us a key control parameter to probe events with 
broader range of initial conditions providing us an opportunity to peep into events with rarer initial conditions which otherwise 
get masked when analysed by centrality binning alone. We find that the inclusion of spectator binning allows one to vary £2 
and £3 independently. We observe that the standard scaling relation between V2 lei and ^ exhibited by centrality bins is 
strongly broken by the spectator neutron bins. Further, the acoustic scaling relation between In (v„/e„) and transverse system 
size is also broken- the strength of the breaking being sensitive to the binning procedure. The introduction of the spectator 
binning allows us to tune over a wide range viscosity driven effects for events with varying initial states but similar final state 
multiplicity. 
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I. INTRODUCTION 


Of all the stages in a heavy-ion collision (HIC), the initial stage is the least understood. However, in order to 
perform a sensitive test of the theoretical framework, e.g. relativistic viscous hydrodynamics, that correctly describes 
the evolution of the strongly interacting matter produced in HIC experiments and therof allows unambiguous extraction 
of the medium properties, e.g. values of the transport coefficients, require a precise knowledge of the initial state (IS). 
It has been shown that depending on the choice of the initial condition that one chooses to evolve the relativistic 
hydrodynamic equations, the value of the extracted shear viscosity to entropy density ratio at the RHIC 200 GeV can 
vary by a factor of 2 [l]-Q . 

The nuclei used in HIC experiments are extended objects. This results in event by event (E/E) geometrical fluc¬ 
tuations in addition to the intrinsic quantum fluctuations of the nuclear wave function. The geometry of the nucleus 
ensures that various characteristics of the IS in HICs like the number of wounded nucleons -/Vp art , number of binary 
collisions IV co n, shape of the overlap region, say the ellipticity £ 2 , are all correlated with the impact parameter b. How¬ 
ever, none of the above IS collision attributes are directly observed in experiments. This makes the job to constrain 
the IS very challenging. The standard method uses the final state (FS) charged particle multiplicity to characterize the 
events into different centrality classes corresponding to different ISs. However, the geometric and quantum E/E fluc¬ 
tuations in the IS result in appreciable variation of b, N pa r t, -/V co n, £2 etc., even within the same centrality bin. Thus, a 
lack of proper knowledge of the IS is a major hindrance towards carrying out precise comparisons between theory and 
experiments. In this work, we focus on the spectators (those nulceons which do not participate in the collision) and 
show that it is possible to extract vital information of the IS by analysing them. The significant role played by spec¬ 
tator asymmetry in the various experimental observables and the possibility of selecting special initial configurations 
in HIC using deformed U nuclei has been pointed out recently (3, Q. In a study based on Monte Carlo Glauber model 
simulations it was suggested that spectator asymmetry could be used to trigger specific collision configurations called 
Body-Tip with sufficient magnetic field and much lower ellipticity which can lead to the disentanglement of chiral 
magnetic effect from its dominant background anisotropic flow in U + U collisions [zj. In Ref. [g] it was demonstrated 
based on A Multi Phase Trasport (AMPT) model simulations that spectator asymmetry could be utilised to identify 
Body-Tip collision configurations in U + U with large forward-backward asymmetry in particle production as well as 
considerably smaller elliptic flow, V 2 ■ In this work, we look into the possibility of probing HICs of non-deformed nuclei 
using spectators. We will in particular focus on Pb + Pb collisions at ^/snn = 2760 GeV. A recent study has pointed 
out the correlation between forward-backward asymmetry in particle production and spectator asymmetry in Pb + Pb 
collisions at ^/snn = 2760 GeV |>j. 

We have constructed an IS observable, namely the total spectator neutron (L+R), which is the sum of the left going 
(L) and right going (R) spectator neutrons that are detected by the zero degree calorimeters (ZDC). The spectator 
protons never reach the ZDC as they are bent by the magnetic field and hit the beam pipe wall in experiments say 
for example at RHIC. Hence we do not consider them. With the present design of the ZDC, faithful measurements of 
spectator neutrons are limited to central and mid-central collisions upto ~ 40% centrality due to clustering effects 0 - 
However, there are suggestions for advanced designs with much improved ZDC performance which could be imple¬ 
mented in the future 11]. In the absence of the IS fluctuations, L + R and b would have a one to one correspondence. 
Both quantities are zero for full overlap collisions and increase for peripheral collisions. The advantage of L + R is 
that it is an observable measured in experiments while b is never measured. We demonstrate that by performing a 
further binning over L + R in addition to the standard centrality binning, it is possible to probe the fireball with novel 
IS conditions as compared to centrality binning alone. This allows us to perform a more accurate comparison between 
model predictions and data, thus enabling a more precise modelling and comprehensive understanding of the response 
of the strongly interacting matter to the IS conditions created in HIC experiments. The outline of the paper is as 
follows: in the next section El we discuss the details of the AMPT model and our simulation as well as our proposed 
methodology to bin events. In Sec. eh we show the main results obtained with this new binning procedure and finally 
summarise in Sec. EO 


II. MODEL 

We have simulated Pb + Pb collisions at ^/snn = 2760 GeV using A Multi Phase Transport (AMPT) 0,10 model 
in the Default version. The AMPT model uses the same initial conditions as HIJING 0 ]. Zhangs parton cascade 
follows to take into account partonic interactions ITH which finally recombine with their parent strings that fragment 
into hadrons within the Lund String Fragmentation model fl6 j] . There is a final stage hadronic afterburner before the 
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FIG. 1. (Color online) Left: The multiplicity distribution for min bias Pb-Pb events at yj snn = 2760 GeV. In alternate white 
and grey bands, different centrality bins are also shown. Middle: The total spectator neutron number L + R distribution for 
the (10 — 15) % centrality. The different L + R bins are also shown. Right: The impact parameter b for the different different 
centrality and L + R bins are shown. 


hadrons freezeout. In this study we have analysed ~ 2 x 10 6 events. 

The standard practice is to first categorize events into different centrality classes according to charged particle 
multiplicity. This is shown in Fig. [I] (a). The different centrality classes are shown in alternate white and grey bands. 
Here we propose to do a second round of binning with the observable L + R within each centrality bin. The L + R 
binning is illustrated in Fig. [T| (b) where the L + R distribution is shown for (10 — 15)% centrality. The L + R 
distribution shows a prominent peak around 85-95 and falls off rapidly on either side - the number of events drop by 
a factor of 5 as L + R shifts by ~ 20. Thus when the analysis is performed based on centrality binning alone, we 
mainly study properties of the events with total spectator neutrons around 85-95. Here with the introduction of the 
L + R binning, we can investigate properties of the rare events with fewer or higher values of L + R compared to the 
centrality mean value. This is the basic reason why additional L + R binning on top of the centrality binning allows us 
to study new IS conditions in HICs. In Fig. |T| (c) we have shown the variation of b with N par t.. We find the centrality 
and L + R bins to follow the same trend. However as shown in Figs. [2j [3] and 01 other IS attributes show different 
trends along centrality and L + R bins that finally translate into different FS obsevable trends. 


III. RESULTS AND DISCUSSION 


We will now present the results of our analysis based on the spectator binning on top of the multiplicity binning. 
We show the results from centrality binning alone by the open traingles joined by a dotted line. Here we propose 
to further bin each such centrality bin into different L + R spectator bins as well. These are shown by the colored 
symbols, each color representing a centrality bin. 

Due to E/E fluctuation of the participant positions, the principal axes of inertia of the participant (P) nucleons 
is shifted as well as tilted with respect to those of the nucleus-nucleus (N) system 03 • Hence, we first perform the 
necessary translation, 

x' = x N -(x N ), ( 1 ) 

y' = VN ~ (vn), (2) 


where (xjvj Vn) denote the nucleon coordinates in the N coordinate system, and further rotate the primed coordinate 
system by the second order participant plane angle ’I / 2 P i so as to coincide with the P coordinate system. is 

obtained as follows, 




ivz. 


E„ 


_/« in# 


Ed 


(3) 


where (r', (/>') is the new 2-D polar coordinate of the ith participant in the primed coordinate system. is obtained 
for n = 2. The initial spatial geometry of the overlap region is encoded in the eccentricities e n defined in Eq. [3] of 
which £2 and £3 will be discussed in this work. 

We will first analyse the distribution of the participants in the plane transverse to the collision axis, which we call 
the collision plane. In Figs. [2] (a) and [2] (b) we have plotted the bin average of the standard deviation in the xp 
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FIG. 2. (Color online) Different qualitative dependence of o^, 
bins. 


(7 y, N co ii and dN c h /drj on iV par t with centrality bins and L + R 


[a x = -^/ (xp) — (xp) 2 ) and j/p ( a y = (t/p) — (yp) 2 ) coordinates of the participants measured with respect to the P 

coordinate system. o x and a y gives us an idea of the initial size of the fireball on the collision plane at the time of 
collision. The curves for cr x and a y vs iV par t show different correlations along centrality and L + R bins. They both 
decrease more rapidly along L + R bins than along centrality bins. This brings us to Fig. [2] (c) where the number of 
binary collisions, N co n is shown normalised to -/Vp ar t, the number of participants. Within a geometry based approach 
of particle production and initial conditions in HICs, e.g. the two component Glauber model, Al part and iV C oii are the 
two most essential ingredients that determine the IS as well as the FS multiplicity [l8l - 5o| . A collision between two 
nucleons (one from each of the colliding nucleus A and B) with coordinates ( xa,Va ) and ( Xb,Vb) on the collision 
plane occurs if they satisfy the following simple geometrical criteria 

{x A - x B ) 2 + (va - Vb) 2 < —(4) 

7 r 

where a^N is the nucleon-nucleon cross section. IV part is the sum of all the nucleons that satisfy Eq. [|] iV C oii 1 on 
the other hand, is the sum of all such possible binary collisions between the participants. The distribution of the 
participants on the collision plane determine the IS eccentricities given by Eq. [3] As shown in Fig. [2] (c), the L + R 
bins of a higher centrality bin mostly have larger value of IV C oii compared to a lower centrality bin. This is so because 
as seen in Figs. [2] (a) and (b), bins with same IV par t but higher centrality occupy a smaller area on the collision plane 
as implied by smaller values of a x and a y . This means the participants are more lined up along the beam axis than 
perpendicular to it, hence having a higher value for iV C oii and smaller a x and a y . Finally in Fig. [2] (d) we have plotted 
dN c h/drj vs N part . Clearly, the correlations along centrality bins is different from that along the L + R bins. In a given 
centrality bin as L + R increases, N part decreases. However, dN^/dr/ almost remains constant which in turn implies 
constancy of the initial energy deposited even though N part decreases. This suggests that in a given centrality, as we 
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go towards bins with larger L + R, the energy deposit pattern changes- one expects to find larger energy gradients 
and more number of energy hot spots (as the same energy is deposited over a smaller transverse area by lesser N pai t). 
This should result in very different viscous effects as one scans over bins with varying L + R at a given centrality which 
will be best borne out by observables on anisotropic flow. 
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FIG. 3. (Color online) Different qualitative dependence of £ 2 , £ 3 , Vi and V3 on N pa ,rt with centrality bins and L + R bins are 
shown. We also show the correlation between £2 and £3 as well as between V2 and V3 with centrality and L + R bins. 


The collective hydrodynamic response converts the IS spatial anisotropy of the fireball as reflected by e n into the 
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FS azimuthal anisotropy of the produced charged particle characterised by the flow observables (v n . V F^' P ), 
dN 

oc 1 + ; 

n= 1 n— 1 

M 


d4> 


pe l + 2 v n cos (n {<j> - ^n P )) = 1 + 2 X! ( Vn ’ x cos ( m + Vn ’V sin ( n ^)) • 


M 


Qn,x = '^2 C0S ( n< &) > Q n ,V = ^2 sin ( n< &) 


Vn,x — ]^j-Qn,xi 


V n,y ~ j\jQn,y, 


Vn = \hl tX +vl y 


( 5 ) 

( 6 ) 
( 7 ) 


Here <j>i is the azimuthal angle of the ith particle, M is the total number of particles and 'F pp is the event plane angle, 
measured using the produced particles 21]. There have been attempts to study the influence of the initial event shape 
on the ensuing fireball dynamics and hence on various FS observables [ 22 l l23j . We now focus our attention on the 
event shape- how L + R binning introduces a control parameter to tune the IS geometry. 
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FIG. 4. (Color online) The correlation between £2 — N co u/N var t and £3 — N co ii/N P art for different centrality and L + R bins. 


From Figs. [3] (a) and (c) it is clear that £2 and £3 can be tuned by simply triggering on different L + R bins within 
a particular centrality bin. Thus the spectator bins allow a direct access to the initial event shape. £2 and £3 show 
almost similar variation along centrality and spectator bins for central events. For mid-central to peripheral bins, 
starting from (20 — 25) % centrality bin, the correlation between £ 2 and fV part along spectator bins is slightly steeper 
compared to that along centrality bins. On the other hand, £3 correlation with N p3iIt is gentler along L + R bins. 
This difference in variation of £2 and £3 along L + R vs centrality bins becomes more transparent in the FS flow 
observables, t >2 and V 3 as shown in Figs. [3] (b) and (d). We note that the L + R bins preserve the ususal linear relation 
between (£ 2 , £ 3 ) and (^ 2 , 113 ). This also gives rise to the £2 — £3 and V 2 — V 3 correlation plots as shown in Figs. [3] (e) 
and (f). It is clear that the spectator binning allows us to access novel geometries in terms of (£ 2 , £ 3 ) pairs which 
are never accessible in centrality binning alone. Moreover, since the trends are quite different along spectator bins, 
they cannot be accessed even if one performs a narrower centrality binning. For (20 — 25) % and more peripheral 
bins, there is much smaller correlation along L + R bins as compared to that along centrality bins, e.g. in case of 
the (20 — 25) % centrality bin while V 3 changes by only 10% for different L + R bins, V 2 varies by around 130%. This 
will allow us to disentangle the effects of V 2 from V 3 on other observables, e.g. the non-linear mode couplings during 
the hydrodynamic expansion that result in correlations in {112 — V 4 ) and ( 112,113 — 115 ) could be studied. In this regard, 
it will be interesting to look at correlation plots of ( 112 , 114 ), ( 112 , 115 ) and ( 113 , 115 ) in data with a combined binning in 
dN c h/drj and L + R bins. Thus, a combined binning in terms of centrality and L + R allows us to disentangle the 
contribution of different IS characteristics on the FS observables. Here we have discussed two such cases, (a) A r part 
and Won and their contribution towards FS dN c h/drj and (b) £2 and £3 and their contribution towards 112 and 113 
respectively. In Fig. |Jwe show the correlation between N co u/N part with £2 and £ 3 . Thus with the introduction of 
L + R bins, we can now study the evolution of similar initial geometry (£ 2 , £ 3 ) but with different mechanism of energy 
deposition (N co u/N part ). 
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Thus studies with L + R bins can complement ongoing studies with <72 bins which also aim at studying the influence 
of the event shape on various FS observables [ 22 l . I23I ] where qi is obtained from the the second harmonic flow vector 

Qi mm 

Qz = \]Ql,x + Q\,y 92 = ( 8 ) 

Here Q 2 , x and Q 2 ,y are computed using particles produced in the FS. It is important to note that unlike the q 2 binning 
procedure, in this method of studying the IS geometry through L + R bins the linearity relationship between (v 2 ,vz) 
and (£ 2 ,£ 3 ) is not essential as the spectators, being IS observables, provide a direct access to the IS geometry. 



FIG. 5. (Color online) vs (1/S) dN c h/dr] for different centrality and L + R bins. The L + R bins break the scaling relation 

between vi/ei and (1 /S)dN c h/dr/ that is exhibited by centrality bins. 



FIG. 6 . (Color online) Acoustic scaling of the hydrodynamic response ln(u„/e n ) vs 1/At with n = 2 (left) and n = 3 (right) 
for different centrality and L + R bins. 


Ideal fluid hydrodynamics is scale invariant and as a result the final value of the ratio (V 2 /£ 2 ) attained is independent 
of the system size [24] . Viscous corrections in a non-ideal fluid arising due to incomplete thermalization introduce 
system size depedence and tend to reduce this ratio. The ratio of the microscopic mean free path A to the typical 
macroscopic system size A, A/A is called the Knudsen number K. /\ -1 is expected to be a good measure of the 
thermalization achieved ,T, 24|. Since develops over some time during which the fireball rapidly expands it is not 
exactly clear what should be the value of A and A to be used to determine the degree of thermalization. In earlier 
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studies, it has been suggested that vi dominantly develops in the early stage of the fireball evolution. Therefore 
A' -1 estimated at time r ~ A t/c s at the onset of transverse expansion of the fireball, should be used as a measure of 
thermalization 0,123- Here is the transversize size of the fireball on the collision plane and c s is the speed of sound. 
It was further argued that A'” 1 at At/c s approximately scales with (1/ S) dN c h / dr) where S is the initial transverse 
area on the collision plane. Similar scaling relation between V 2 /E 2 and (1 /S)dN c h/dr] is also expected in the low 
density regime [25F27]. Previous analysis of experimental data with centrality binning have indeed found very good 
scaling relation between V 2 /E 2 vs {l/S)dN c h/drj for different systems like Cu+Cu and Au+Au Hi- The detailed 
mechanism of energy deposition in the intial stages of a HIC event which ultimately leads to particle production is 
yet to be understood completely. With centrality as the only tuning parameter to separate different initial conditions, 
it is difficult to discriminate between predictions from models with different mechanisms of particle production. Now, 
with the introduction of L + R bins, we are able to pin down the initial conditions more precisely. In Fig. [5] we have 
plotted the ratio ^ 2/^2 vs (l/S)dN c h/dr) with centrality and L + R bins. The centrality averaged points exhibit the 
usual trend of an inital fast rise and final saturation of the ratio V 2 /S 2 with (1 /S) dN c h/dy. However, the L + R bins 
in each centrality show the opposite behaviour and breaks the usual scaling relation. In a given centrality bin, with 
rise in L + R, the transverse overlap area S sharply falls while dNch/dy is almost constant. Thus (1 /S) dN c h/dy which 
acts like a proxy for K _1 increases with L + R although the hydrodynamic response V 2 IS 2 falls sharply. As mentioned 
earlier, we expect more (less) hot spots and gradients in the initial energy profile of events with larger (smaller) L + R 
leading to more (less) viscous correction. This ultimately leads to inefficient conversion of the initial e n to the final 
v n in events with larger L + R as compared to those with smaller L + R. 

The ratio v„ /e„ is also expected to exhibit acoustic scaling and receives viscous corrections that grow exponentially 
as n 2 and 1/A T [H® 


In 



4 n 2 ?? 

3 A t Ts ’ 


( 9 ) 


with the typical initial transverse size of the system A^ given by 1 /Aj- = yC/cr^ + 1 /(jy. Such acoustic scaling was 
found in data across a wide range of beam energies and the shear viscosity to entropy density ratio, 77 /s was extracted 
from the slope of the plot, In vs 1/Aj- [29l - l3T| . These studies were done with centrality binning alone. Here we 

perform a consistency check of such scaling laws using simulated data by including L + R bins in every centrality bin. 
In Figs. ED (a) and (b) we have plotted In ( 772 /^ 2 ) and In ( 773 / 53 ) respectively vs 1/A-r- We find approximate simulated 
data collapse for (0 — 40) % centrality as well as their corresponding L + R bins. On a closer scrutiniy, it seems that 
the slope parameter for bins of different L + R but a particular centrality is different from those of a different centrality 
resulting in a mild breaking of the acoustic scaling. Thus, the introduction of the L + R bins enable a more refined 
extraction of y/s from data. 



FIG. 7. (Color online) Different qualitative dependence of a x , N co u and dN c h/dri on A part with the reverse binning procedure 
(first binned by L + R followed by dN c h/drj) are shown. 


So far we have discussed results for events first binned by dN c h/dy followed by L + R. If b was the only E/E 
fluctuating quantity, then 6 , dN c h/dy and L + R would have a one to one correspondence and hence the final results 
would be independent of the order of the binning procedure. However, as discussed earlier, in HICs there are additional 
sources of E/E fluctuations apart from the geometrical fluctuation in b. This means that the final results are sensitive 
to the order of the binning procedure. In order to illustrate this point we have also analysed the events in the reverse 









FIG. 8. (Color online) vi/si vs (1/S) dN c h/dri and 1 /At with the reverse binning procedure (first binned by L + R followed by 
dNch/dr/) are shown. 


binning procedure: first we bin by L + R followed by dN c h/dy. We have shown a few results in Figs. [7] and [5] As seen 
from Fig. [7] (a), cr x which is also a measure of the initial system size is almost constant in a particular L + R bin even 
though N part changes resulting in the strong variation of N co u/N part within a L + R bin as shown in Fig. [7] (b). This 
finally translates into a stronger variation in dN c h/drj within a L + R bin with N part as compared to the variation 
between different L + R bins. This further results in the trends between v-i/ei and l/SdN c h/dy as seen in Fig. [8] (a). 
Within a particular L + R bin, the bins with higher l/SdN c h/dy have higher N co u/N part ratio than the average trend 
resulting in smaller 112/62 compared to the average trend. This trend was seen even in the earlier binning procedure. 
As seen in Fig. [5] (b), the acoustic scaling relation as in Eqn. [9] does not hold anymore as the slope for the trend of 
the average L + R bins is much softer than the slope along different dN c h/dy bins in the same L + R bin. Thus this 
breaking of the acoustic scaling relation which was mild in the earlier binning procedure becomes much stronger in 
the reverse binning proccedure. 

So far all our results have been from the AMPT model. We will now show a few results for the combined binning 
procedure with dN c h/dy followed by L + R in HIJING event generator. In Figs. [9] (a) and (b) we show the results 
for b and dN c h/dy and find them to be very similar to what we obtain in the case of AMPT in Figs. [U (c) and[U] (d). 
AMPT takes into account later stage interactions in the partonic as well as hadronic phases while in HIJING such 
interactions are absent. As seen in Fig. [lUl this results in very different trends for (pt) with N part in AMPT compared 
to HIJING. In AMPT, we find the (j>t) trends along L + R bins to be much different to that of the average trend 
along different centrality bins with ~ 10 — 15% variation in a particular dN c h/dr] bin for different L + R bins. On the 
other hand, in HIJING the (pt) trends do not depend on the binning procedure and is a single-valued with N part . 
We conclude that this starkly different trends of (pr) in the two cases should be stemming from the fact that medium 
effects which are taken into account in AMPT are missing in HIJING. In the different L + R bins there is different 
degree of medium interactions and collectivity resulting in different values of the final (pt)- Our earlier observation of 
different values of V2/62 along different L + R bins support the above conclusion as well. Thus measurement of (j>t) 
with such combined binning procedure can also throw light on the degree of collectivity achieved. We note from Fig. [5] 
that in a particular dN c h/drj bin, those L + R bins which have the least 1)2/62 also have the highest (pr) in Fig- EH 
(a). 


IV. SUMMARY 

We have demonstrated using the AMPT model the important role played by the spectators to determine the initial 
condition in heavy-ion collisions. The standard procedure involves binning events by their final state multiplicity. 
This however puts events with varying inital conditions into the same bin as long as they produce similar multiplicity. 
We demonstrate that by further binning events according to the total number of spectator neutrons, it is possible to 
separate events with different initial conditions which were earlier clubbed together in the same centrality bin. This 
new methodology provides an opportunity to study events with rare initial conditions. Moreover it is possible to vary 
62 and £3 independently of each other. This enables one to extract the contribution due to non-linear mode coupling 
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FIG. 9. (Color online) The trends for dN c h/dr/ and b with Apart as obtained in HIJING where the events were binned by 
dN c h/drj followed by L + R. 
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FIG. 10. (Color online) The trend for (pr) with Apart as obtained in AMPT (left) and HIJING (right) where the events were 
binned by dN c h/dr/ followed by L + R. 


between V2 and V4 and (112,113) and V5. It is important to note that for this purpose it is not essential to know L + R 
very accurately. We found the variation of dN c h/dri with Af part to be much different for L + R bins compared to usual 
centrality bins thus allowing us to study different energy deposition mechanism within the same centrality. We argued 
that in a given centrality bin, larger L + R bins have higher energy gradients and more number of energy hot spots 
as compared to smaller L + R bins which result in strong inhomogenities in the inital conditions. This calls for larger 
viscosity driven effects and hence smaller v? /£2 for bins with higher L + R. This also results in the breaking of the 
scaling relation between V2/S2 and j,dN c h/dr] by the L + R bins. A comparatively milder breaking of the acoustic 
scaling relation between In (v n /e n ) and the initial system transverse size is observed for both centrality as well as 
L + R bins. The results from this study suggest that one might be able to extract a more accurate value of the 77/s 
ratio with the introduction of the L + R bins. We also observe that (j>t) in combined bins of dN c h/di~i and L + R is a 
good probe to measure the degree of medium interaction and again in this case a precise measurement of L + R is not 
necessary. Hence even with the current performance of the ZDCs, we should be able to perform some of these analysis 
in data. 
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